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' The isoscalar and isovector collective multipole excitations in stable nuclei are studied in the 

o ■ 

' framework of relativistic random-phase approximation with the vacuum polarization arising from 

(N 

■ the nucleon-antinucleon field. A fully self-consistent calculation which guarantees the decoupling of 

the spurious state and the conservation of the multipole-transition current is carried out by using the 

' derivative-expansion method for the description of the vacuum contribution. A remarkable effect 

^S) ' of the inclusion of the vacuum polarization is the increase of the effective mass, m*/mjv ~ 0.8; 

> 

for all multipole modes, the energy-weighted sum rule values with the vacuum polarization are 
' smaller than those of the relativistic model without the vacuum polarization, which typically has 



o 
o 



the effective mass m* /niM ~ 0.6. Also, the present model can give an excellent agreement with 
experimental data on the excitation energy, in particular, for the isoscalar quadrupole resonances 
in which it was previously reported that the calculated energies in the relativistic model are about 
1-2 MeV above the experimental values. It is shown, further, that the change of the shell structure 
due to the inclusion of the vacuum polarization plays an important role in the improvement of the 
discrepancies seen in the dipole compression modes. On the other hand, the isoscalar monopole 



X 

' resonance has a similar peak whether or not the vacuum polarization is taken into account, if the 

compression modulus is kept the same in the analysis. 
PACS immber(s): 21.10.-k, 21.60. Jz, 13.75.Cs, 23.20.-g 

I. INTRODUCTION 



The giant resonance of nuclei is of great interest because its energy can be directly related to the 



compressibility of nuclear matter 



which places important constraints on the theoretical description 
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of nuclear properties, supernova explosions, neutron stars, heavy ion collisions, and so on. The reliability 
of the model employed for the giant-resonance analyses may be thus confirmed by the comparison with 
the experimental energy-weighted sum rule (EWSR) distributions. On the other hand, only the ratio 
of the EWSR between different nuclei, rather than the EWSR strength, was available because the sum 



Q 



rule strength was dependent on model assumptions, as pointed out in Ref. [3j. A knowledge of the 
magnitude of the EWSR for the electromagnetic response in addition to its distribution, however, will 
impose a further constraint with regard to the nuclear model, and offer valuable information for the study 
of a two-photon exchange process through the nuclear particle-hole excitations such as the dispersion 
correction to the elastic electron scattering 4J and nuclear-polarization correction to the muonic levels j^, 
which are approximately proportional to the EWSR strength. 

The comparison between the experimental data and the theoretical analysis on the nuclear excitations 
provides a critical test of the quantum hadrodynamics (QHD) in a many-body systemy|. A standard 
tool for the theoretical description of nuclear excitations is the random-phase approximation (RPA). In 



particular, many studies have been done by u sing relativistic RPAlKKJ 

2, 3 3 3 u, 3 3 M mH 17, 3 19, 20, 3 



field (RMF) basis 



(RRPA) with the relativistic mean- 



1. On the basis of the works of 



the last decade, it has been realized that the self-consistent method with nonlinear effective Lagrangians 
is essential for a quantitative description of excited states. In this context, it is very interesting that 
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the negative-energy nucleon state in the mean-field basis seriously contributes to excitation energies 
the current conservation of the transition densities 3 3 3 > ^"^^ ^^'^ decoupling of the center-of- 
mass motion 

1313 

Furthermore, it was found that the negative-energy RRPA eigenstates, generated 
from the RRPA equation with the fully consistent basis including negative-e nerg y RMF states, have a 
significant role in gauge invariance of the response to an electromagnetic field |22l|. 

The success of the relativistic nuclear model is mainly ascribed to the small effective mass, m* /ttin ~ 
0.6, which yields the appropriate single-particle splittings in a nuclear ground state. The fact that the 
centroid (or peak) energies and the distributions of the nuclear excitations in the isoscalar monopole mode 
and the isovector dipole mode are well reproduced [loj seems to support the validity of the small effective 
mass. However, the physical quantities which are more sensitive to the variety of the nucleon mass are 
the strength of the response function and the centroid energy of the giant-resonance mode which the 
shell structure around the Fermi surface affects significantly. From this viewpoint, an interesting fact in 
the recent study of Ref. 23] is that the /?-decay lifetimes require an enhanced effective mass to reproduce 
the data. 



The small effective mass is caused by ignoring negative-energy nucleons (or antinucleons). If the Dirac 
sea filled by negative-energy nucleons is duly taken into account, then the effective mass converges at an 
enhanced value naturally because the Lorentz scalar field is strongly suppressed due to the feedback effect 
from antinuclcon states(2^. Inclusion of the Dirac sea, however, raises a problem in fitting the single- 
particle energies, requiring a new mechanism to supply the spin-orbit force. A simplest solution would be 
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In this work 



to add the tensor-coupling term for the meson- nuclcon interaction to the Lagrangian 
we employ a new extended relativistic model with tensor coupling for the vector meson incorporating 
the Dirac sea effects of nucleons. We emphasize that the self-consistency is kept throughout the RHA + 
RRPA calculation, as verified by obser ving the spurious state in the isoscalar dipole mode and the current 
conservation of the transition densities 23. We present a systematic study of the giant monopolc, dipole, 
and quadrupole resonances for the different effective mass in the fully consistent RRPA framework. 

This paper is organized as follows. In Sec. II we introduce the effective Lagrangian with renormalized 
negative-energy nucleon contribution by using the derivative-expansion method. With this approach, the 
contribution from negative-energy nucleons to the effective action can be expressed by the functional of 
meson fields. Such a action gives the meson propagators dressed by the vacuum polarization of nuclcon- 
antinucleon pairs automatically. The details of the computational procedure for the self-consistent RHA 
+ RRPA calculation with the vacuum correction are described in Sec. III. The vacuum effect to the 
decoupling of the center-of-mass motion and the way to conserve the RPA transition densities in the 
model, including the tensor-coupling interaction, are also shown here. In Sec. IV, the numerical results 
of the present calculation for the Coulomb response function in the spherical nuclei are compared with 
those of the TMl and NL3 parameter sets. In particular, the giant monopolc, dipole, and quadrupole 
resonances of the ^°Zr, ^^®Sn, and ^°*Pb as well as their energy-weighted sum rules are discussed in 
detail. Finally, we give a summary of our calculations and an outlook of future applications in Sec. V. 

II. FORMALISM 

Within a relativistic framework, it is inevitable that the negative-energy nucleon states contribute 
to the nuclear structure as a vacuum fiuctuation. Recent investigations of the RRPA in the QHD 
model, however, have neglected the actual antinucleon degrees of freedom; the basis set used for the 
RRPA calculation is prepared by the RMF theory wherein only positive-energy nucleons are taken into 
account, so that the nuclear-polarization tensor in the RRPA equation is described without the vacuum 



polarization of the nuclcon-antinuclcon field. On the other hand, the rigorous calculation of basis set 
under the full one- nucleon- loop contribution, referred to as the relativistic Hartree approach (RHA), 
has been developed by using the Green function method ^ 2^ 2^. With this technique, however, the 
estimation of the vacuum-polarization effect is too hard even for the ground state of the light nucleus 
to extend it to the RRPA formalism including the particle-hole correlations, especially for heavy nuclei. 



Fortunately, it was found in Ref. 



that the lowest-order contribution of the derivative expansion 



from the one-baryon-loop correction to the effective action is in excellent agreement with the rigorous 
result. The derivative-expansion method describes the vacuum contribution by the additional new terms 
of meson fields in the effective Lagrangianjs^. As a result, it makes it easy to perform the RRPA 
calculation with vacuum polarization because the meson propagators are constructed by including the 
vacuum part of the nuclear response as well as self- interaction terms 27 1 . 

A remaining problem in nuclear models with the vacuum-polarization effect was that the magnitude 
of the spin-orbit splittings in the shell structure of the mean field is insufficient in comparison with 
the observed data because the vacuum-polarization effect decreases the spin-orbit potential of nucleons. 
Reasonable results for the ground-state properties of the binding energy and the charge radius as well as 
the spin-orbit splittings were given by Mao 2g , who has investigated the single-particle levels of nucleons 
and antinucleons by incorporating the effect of the tensor-coupling meson-nucleon interaction. As far as 
the spin-orbit splittings are concerned, the tensor coupling of meson-nucleon interaction, which directly 
generates the spin-orbit force in the RHA calculation, is considerably effective for the simple extension of 

n 

the model in comparison with other candidates such as the surface pion condensation 31] and the Fock 
contribution of the pion exchange. Hence, in the present paper we shall also incorporate the tensor- 
coupling interaction for the w meson. Although the model is now non-renormalizable and the divergence 
should be removed by taking an infinite number of couplings into account, we assume that such a vacuum 
contribution from the tensor-coupling interaction is negligible. The following Lagrangian density is thus 



employed for the present model; 
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(1) 



where ip is the Dhac spinor of the nucleon, and cr, w^, p^, and represent the isoscalar-scalar meson, 
isoscalar-vector meson, isovector-vector meson fields, and the electromagnetic field, respectively. Here 
<^tiu = (*/2)[7p, 7iy] is the spin tensor, Q = (1 + T3)/2 projects the proton, and Ua{cr) — ■^g2<^^ + jids'^'^ 
{Uu,{iLJf_i) ~ jc^i^fi^^)'^) denotes the self interaction of isoscalar-scalar (-vector) meson. The superscript 
+ in the nuclcon-field operators in Eq. |^ means that only the positive-energy states are to be treated 
explicitly. The functions Vpicr), Zp{a), Z %(a), a nd Zpia) which describe the vacuum effect of nucleons 



are obtained by well-known prescriptions 
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The paths of classical meson and photon fields are determined under the condition that the contribution 
from the first order of the quantum fiuctuations vanishes in the action. Assuming the stationary and 
spherical system, the classical fields ct, (Dq, po.3, and Aq satisfy the following coupled equations: 
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V'Ao - e(^+Q7>+) - V • (Z#(a)VAo), 



(6) 
(7) 
(8) 
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The contributions from negative-energy nucleon states to source terms are contained in Vp and Zp- In 
a mean- field approximation, nucleon states satisfy a Dirac equation under the classical fields a, ojo, po,3, 
and Aq: 
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Tpa = l^alpai (10) 



for which there are solutions with negative energies as well as ordinary positive energies. The RHA basis 
consists of the positive- and negative-energy eigenstates obtained by solving the Dirac equation ()10|) and 
the Laplace equations ©-0 iteratively, until convergence is reached. 

Since the vacuum-polarization effect is taken into account in the RHA basis, the RPA diagram also has 
to involve the vacuum contribution to retain the consistency. This is achieved by solving the following 
Bete-Salpeter equation with the full- unperturbed polarization tensor Hh- 

^flPA(^^ r^p, q- uj) =nH(r^ r^p, q;Lu) 

+ E / ^-^nHiT\r-,p,k,;co)D^{k^,k2;co)nnPAiT\T'';k,,q;iu), 

(11) 

where the summation is over the meson fields, and denotes the meson propagators with coupling 
constants: 

i?0(fci,fe2;c^) = - ^ ^i2nf5{ki - k^) 
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UJ^-kl-mlJ (27r) 



/ -^U'Jia){k^,k,)D'^{k,,kr,u;), (12) 



Dl{ki,k,; 0.) ^-^-^—^i2nmk, - k,) 

D%k^M;^)^^^ ^(27r)3j(fci-fc2), (14) 

D\{kiM;^) = , 2 (2^)^(5(fci ^fca), (15) 



where the cr- and w-meson propagators have to be constructed with the self- interaction term|15|. T's in 
Eq. are the 4x4 matrices which denote the vertex couplings, 

rel,7p + ^'5''a^„^7M-27M- (16) 



to the cr, a;, p mesons, and photon, respectively. The concrete expression of Hh is given by. 
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UJ — {uJa — ujf}) + iri/2 UJ — {uJa — ujp) — iriji' 
where ei? is the Fermi energy. They are sketched in Fig. where II^i is composed of the ordinary 
particle-hole excitation and the Pauli-principle violating excitation between the positive-energy and the 
negative-energy nuclcon states. On the other hand, Hi? is composed of the excitation between the 
positive-energy particle and the negative-energy hole as the Fermi sea is unoccupied. The latter is 
divergent, and therefore the renormalization procedure is required to extract the physical contribution 
from it. To our knowledge, IIj? has never been calculated explicitly in a finite system and it is still a 
remaining task. However, the use of the derivative-expansion method, which renormalizes the divergence 
in the effective action, enable us to evaluate the vacuum contribution to the polarization tensor within 
this approximation. As indicated in our previous paper 271, the second-order fiuctuation terms of meson 
fields in the effective action give the meson propagators including the vacuum polarization expressed 
by the functions of meson fields. For the RPA equation we shall define the finite unperturbed 

polarization tensor as 
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where the vacuum-polarization part Hj? is given as 
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where the superscript in Hp indicates mesons at initial and final vertices. We have shown in Ref. 
that the leading order of the derivative expansion in RHA calculation reproduces the numerical results of 
the exact Green function method very well. It is expected that also in the RPA calculation, it is sufficient 
to expand up to the leading order of the derivative terms, as far as the low-momentum transfer region 
are concerned. The verification of the self-consistency will provide a check that the vacuum polarizations 
in the RPA level are sufficiently described by Eqs. I|21|) - (|27|) . In the next section, we will show that the 
spurious state in the isoscalar dipole mode is decoupled and the transition density is conserved with 11^? 
of Eqs. (PTl)-(|77|l. 

III. NUMERICAL DETAIL OF VACUUM CONTRIBUTION, TENSOR-COUPLING 
EFFECT, AND ANTINUCLEON SOLUTION 

We shall explain here the main points of the calculation of the self-consistent RHA + RRPA, including 
the vacuum polarization. First of all, we perform the RHA calculation, which gives self-consistent mean- 
field potentials by making the iteration of the nucleon- and meson-field calculation. The parameter 
sets used in the present work are listed in Table I, where the Fermi momentum, the binding energy, 
the effective mass, the compression modulus, and the asymmetry energy in nuclear matter are also 



shown. TMl 



mass, tne compr 
3 and NLslsa ? 



and NL3|37| are the well-known parameter sets which give good descriptions of nuclear 
properties and RHATl, RHAT2, and RHAT3 are the parameter sets employed for the present model 
explained in the previous section. They are determined by firstly surveying the parameters in nuclear 
matter and choosing the sets which have the compression modulus K = 270, 281, and 325 MeV and 
then, by fitting to the total binding energies, charge radii, and single-particle energies of the spherical 
nuclei in the RHA calculation, with a fine adjustment of the coupling constants of the a and w mesons 
in addition to the coupling constants of the p meson and of the w-nucleon tensor coupling. As seen in 
Table II, they reproduce the nuclear ground-state properties as well as TMl and NL3. 

The mean-filed potential of each parameter set is then used to generate the basis set of nucleons. The 
RHA basis functions in the RRPA calculation are obtained by sol ving the single-particle Dirac equation 
(|10|l using the method of discretization for the continuum states[l^. For the excited states, we have 
included the positive-energy states up to 200 MeV and practically all of the negative-energy bound states 
(up to —1875 MeV), which contain a large number of bound states in the negative-energy spectrum, since 
the scalar and vector potentials add in the effective negative-energy potential. The RRPA equation (|ll(l 
is solved directly by the inverse matrix method in the momentum space which allows the inclusion of 



as many configurations as required without changing the size of the matrix. Figure [21 shows the RRPA 
Coulomb response of the isoscalar dipole resonance in 208pj-, -^yj^jj ^j^g parameter set RHATl. There, 
we choose the momentum transfer, q = 10 MeV, and the width, rj — 1 MeV. The sohd curve is the 
complete result and the dashed one is the result without Hp"^ (Eq. (|21|l ). ignoring a part of the vacuum 
polarization, in the RRPA calculation. We see that the position of the strength peak which is composed of 
the ccnter-of-mass motion gets shifted to zero energy, as expected, in the fully consistent RHA + RRPA 
calculation. For the dashed line, on the other hand, we find that the spurious component is removed 
insufficiently in the real excitation-energy range. Thus, we conclude that the inclusion of the vacuum 
polarization defined by Eqs. (|21|) - (|27|l in the RRPA equation is important for the self-consistency. 

Another important aspect for establishing the correctness of the self-consistent RRPA calculations is 
to check whether transition charge density p\ {q) and current density jxL (q) connectin g th e ground state 



and the excited states for different multipolarity A and L, satisfy the conservation law 
lun denoting the excitation energy of the nucleus, the conservation relation is given by 



g th e grounc 



With 



i^NPxiq) = -q\j ^ hx-i{q)+q\j ^'^^''^ jAA+i(g)- (28) 

The violation of the current conservation in the transition densities, defined by the difference between 
the left-hand and right-hand sides of Eq. H28|) (hereafter referred to as Ap ) should be very small with 
a large configuration space. However, the RRPA states calculated with the unperturbed-polarization 
tensor defined by Eq. H17|) do not satisfy the conservation law as shown by the solid line in Fig. 13 This 
violation is due to the tensor-coupling term for the u meson. In order to make the transition current 
conserved, the additional term is required in the photon-w meson polarization tensor (More detailed 
argument is given in the Appendix A.); 

^^,,(Q7^7'' + ^ (0|r[j^(a:)j-(y)] |0) + (0|[j^;(a:),£^(2/)]o|0), (29) 

where [• • • ]o denotes the equal-time commutation relation, and 

j'X^i^Qi'^i^. (30) 

jl^^^'^^ + J^d^ii^a'-^^k), (31) 

£a = ^r^Via'^V. (32) 
2m AT 



The second term of the right-hand side of Eq. (|29|l is called the "Schwinger term'' 



. In the present 



case, the Schwinger term exists in the transverse polarization tensor. Including this contribution, the Ap 
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shown by the dashed hne m Fig. |3| beeomes neghgibly smaU, that is, the conservation law is achieved. 
Since the complete set of the basis is assumed in the proof of the conservation of the RPA transition 

.. n . 

densities [13|, this result also confirms that our truncation in the spectral method for the positive and 
negative energies is reasonable. 

Another, and main role of the tensor-coupling term in the present nuclear model is to produce the 
proper spin-orbit force in the single-particle basis. For the mean-field calculation, we do not discuss it 
because Mao has already investigated its effect in detail 2^ . Here, we shall briefly mention the tensor- 
coupling effect appearing in the RRPA states. Since the tensor-coupling term restores the spin-orbit 
force with a large effective mass, the RRPA states may be influenced by the change of the single-particle 
shell structure in addition to the residual tensor-coupling interaction. For the collective modes, however, 
the calculations with and without the tensor-coupling interaction give the similar result with each other, 
if the parameter sets producing the similar compressibility and the similar effective mass are used y. A 
large effect due to the tensor-coupling term appears in the low-lying states, in particular, in the higher 
multipole mode as the quadrupole and octupole excitations, in which the shell structure of the Fermi- 
energy region concerns directly. This fact means that the tensor couplings enable us to fit the energy 
spectrum of the low-lying states phenomenologically without changing the giant-resonance structure. 

In the last part of this section, we shall discuss the RRPA cigenstates with the negative energy. As 
discussed in Ref. [2^, the fully-consistent RRPA equation has the solutions for the cigenstates with the 
negative energy which represent the blocking effect of the nucleon-antinuclcon creation due to the states 
occupied by the spectator nucleon. In the present model with the vacuum polarization, a negative-energy 
solution is interpreted as an antinuclcon one and also is composed of a part of the complete set of the 
RRPA states together with the ordinary positive-energy solutions. As an example, the longitudinal 
(Coulomb) and transverse form factors for the isovector lowest-antinucleon dipole state (1507.2 MeV) 
with the parameter set RHATl is depicted in Fig. 21 The interesting thing seen in this figure is that 
the transverse form factor is considerably larger than the Coulomb form factor. It is also seen that the 
form factors of antinuclcon states have a peak in a very low momentum region than the corresponding 
transition energy, because the states are bound strongly. These facts are very similar to those in the 
RRPA without the vacuum polarization except for the position of the excitation energy; the potential 
in the antinuclcon sector becomes extremely deep by ignoring the vacuum polarization, so that the 
lowest-antinucleon state is located at a lower energy than in the present case. 
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IV. RESULTS AND DISCUSSION 

We now discuss the calculated results of the Coulomb response function for collective excitations 
predicted by the relativistic nuclear model with the vacuum polarization. Theoretically, the nucleus 
cannot be constructed without the large effective mass if the vacuum of nucleons is duly accounted 

nn 

for[2^|40|. Although the large effective mass changes the shell structure obtained by the models without 
the vacuum polarization, the nuclear ground-state properties come out to be similar for all the three 
parameter sets used in this work. For nuclear excitations, however, we can see the vacuum-polarization 
effect through the effective mass effect in the EWSR strength, and in the excitation energies of some 
giant resonances. 

A. Energy- weighted sum rules 

Energy-weighted sum rules depend on the completeness of the particle-hole basis and are in some 
case preserved in the RPA. Therefore, the EWSR is often used as an indicator of the strength of the 
nuclear response in the calculation which requires the summation over all nuclear intermediate states. 
The EWSR value for the Coulomb response is given by 

Y^B{E\;Q^l')ojp'=" ^f^[mp^)jujdujlmURPA{Ql'>,Ql'';q,q;co) (A ^ 0), (33) 
/' 

Y,B{EO;0^r)u;r'^° 1^ (^j,^ J u;du;lmURPA{Ql',Ql';q,q;u;) (A = 0), (34) 

where A denotes the multipolarity. In the present calculation we take a small momentum transfer q = 10 
MeV. The results for A = 0, 1, and 2 of ^°Zr, "''Sn, and ^°^Ph are shown in Table III, where we take the 
summation over the RPA states up to 70 MeV for the monopole states and 50 MeV for the other states. 
From this table, we find that all relativistic results of the EWSR value are somewhat larger than the 



classical EWSR values of Ref. 



41 1 in any multipole states. The reason has been previously presumed as 
4^. 



being due to the effective mass jl3il42l |. Certainly, the EWSR values calculated from the RRPA states are 
approximately proportional to the effective mass as depicted in Fig. [3 The model including the vacuum 
polarization cannot shift the effective mass to less than m* /niN — 0.75 due to the feedback effect from 
the vacuum. As a result, the EWSR's are to be suppressed in comparison with the value calculated 
with TMl and NL3. Nevertheless, they are still enhanced by 10 — 30% over the nonrelativistic classical 
results. In fact, the excess strength has been reported in the analyses of the giant monopole resonance, 
where the single folding (deformed potential) model led to the strength with 154 — 160 (112 — 129)% of 
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the EWSR in the angular distribution of the differential cross section for the 13.7 keV peak in ^°*Pby|. 
The corresponding EWSR for the monopole states in ^"'^Pb with RHATl (134%) is close to the analysis 
with the deformed potential model, while that with NL3 (147%) is close to the analysis with the single 
folding model. 

Although the present results in Table III are thus consistent with the scattering analysis, this model 
dependence of the results suggests that the validit y of the nuclear model with vacuum polarization has 
to be verified by other experimental data. In Ref . |23| , /3-decay rates of r-process in neutron rich nuclei 
were systematically investigated and it was found that, in order to explain observed data, both of the 
relativistic effective mass, called the Dirac mass there, and the empirical effective mass derived from the 
non-rclativistic analysis, which is related to the time like component of the vector interaction, should 
be considerably larger than those of the typical relativistic models. With an additional isoscalar tensor- 
coupling interaction, the relativistic effective mass m* /m = 0.67, which is considerably higher than that 
of the typical relativistic models, gave a significant improvement in the /?-decay analysis. Although the 
experiment demands a larger value, however, the effective mass could not be increased more, provided 
that the ground-state properties of finite nuclei, including the spin-orbit splittings, are reproduced at 
the same time. The present nuclear model, in contrast, generates the relativistic effective mass larger 
than m* jm — 0.75 by means of the vacuum-polarization effect. Inclusion of the vacuum polarization, 
hence, may resolve these difficulties and support the validity of the present model. In addition to the 
/9-decay analysis, it would be interesting to perform the analysis of anomaly in the nuclear-polarization 
corrections for the ^"Zr, ^°*Pb, and Sn isotope by the present extended relativistic nuclear models; the 
nuclear-polarization correction is considerably sensitive to the EWSR strength and the anomaly seems 

4^. 



to require the effective mass effect 



22, 



The vacuum effect will be verified more through such a 



systematic study of the effective mass in the relativistic model. 

B. Giant monopole resonances 

The compressibility of nuclear matter is important in the description of the properties of nuclei. It 
is generally accepted that the best procedure to determine the compression modulus is to calculate 
isoscalar monopole energies regarded as the "breathing mode" of the nucleus by using the microscopic 
models. In the present analysis we examine the predictions of the various parameter sets of the effective 
Lagrangian listed in Table I within the framework of the self-consistent RRPA. The isoscalar-giant 
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monopole-resonance (ISGMR) strength distributions for ^"Zr, ^^^Sn, and ^''^Pb are displayed in Fig. 
where we set the momentum transfer q ^ 10 MeV and the width rj — 1 MeV. The sohd curve is the result 
with the parameter set RHATl and the dashed (the dashed-dotted) one is the result with the parameter 

M 3- It 



set TMl (NL3). The experimental data of excitation energies is also indicated by the arrow 
is found that the difference in these figures among the parameter sets is small. This means that the 
inclusion of the vacuum polarization has a negligible effect in the position of the monopolc excitation 
energy if the compressibility is fixed to a certain value. In Fig.[7| we show the ISGMR centroid for these 
nuclei as a function of the compressibility. Evidently we see that the centroid energies with the similar 
compressibility coincide with each other, and are irrelevant to the inclusion of the vacuum polarization. 

Youngblood et al. obtained the compressibility K ~ 220 — 240 MeV from the analysis of experiment 
with the various non-relativistic calculations 44] , while the analysis with ordinary relativistic models 
gave K — 250 — 270 MeV^^. The origin of the difference between the relativistic and non-relativistic 
predictions is not understood so far. The present analysis including the vacuum polarization did not 
resolve this discrepancy as far as the results of 208p|-, ^^^.g concerned. However, the present results for ^''Zr 
and ^^^Sn indicate that the relativistic model may require the lower value for the compression modulus. 
For this problem we need further investigations including the pairing force. 

In contrast to the ISGMR, the studies to locate the isovector giant monopole resonance (IVGMR) have 
led to ambiguous results because it is in general found at higher excitation energy with less collectivity. 
The calculated IVGMR strengths induced by mainly the p-meson interaction in the present models are 
also distributed widely in fragments. We display the IVGMR strength distributions for ^"Zr, ^^^Sn, and 
^ospb in Fig.lSlwith the three parameter sets RHATl, TMl, and NL3. The resuhs do not depend much 
on the particular model used, and t he p ositions of the IVGMR strength are in good agreement with the 
experiments indicated by the arrow |46l l47 |. 

C. Giant dipole resonances 

The isoscalar-giant dipole resonance (ISGDR) contains the spurious component caused by the center- 
of-mass motion. In the previous section, we indicated that the spurious state appears at zero excitation 
energy with the fully self-consistent calculation, but the physical excitations with q ~ 10 MeV become 
negligibly small since the spurious state has an exceedingly strong collectivity at low-momentum transfer. 
In order to identify the physical excitations, hence, we employ a larger momentum transfer, q — 100 
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MeV, in the ISGDR mode. The ISGDR strengths for ^^Zr, "^gn, and ^ospb are depicted in Fig. M 
ExperimentaUy, it is observed to spht into upper- and lower-energy components in the ISGDR strengths 
(the experimental peaks are indicated by two arrows in the figure [45I 4^). Both in the model with 
(the solid curve) and without the vacuum polarization (the dashed and dash-dotted curves), we see that 
there are two components similar to observed ones. In comparison with the experiments, however, the 
splittings between the upper and lower peaks are somewhat large in all calculations employed here. On 
the other hand, we find that calculations reproduce the peaks of the lower-energy region in ^^^Sn and 
^"^Pb. Moreover, the results with the vacuum polarization for the upper component are different from 
other models, and approach to the data, although the predicted peaks are still 2 — 5 MeV higher than 
the experimental ones. It has been reported that the upper component is the compression mode, whose 
energy depends on the compression modulus, while the energy of the lower component changes with mass 
but is essentially independent of compression modulus 2o[ Li^l. However, in the present model with the 
parameter set RHATl, which has a similar compression modulus to that with TMl and NL3 (see Table 
I), the energies in upper region for all nuclei certainly shift to the lower energy than those with TMl and 
NL3. This suggests that not only the compression modulus but also the shell structure is important to 
resolve the discrepancy in the position of the peak for the upper-energy component. 

The isovector giant dipole resonance (IVGDR) are also illustrated in Fig. ^| There, we find that 
due to the effective mass suppressed by the vacuum polarization, the strengths of the RHATl (the solid 
line) become weak for all nuclei. The centroid energy of ^'^Zr with RHATl (16.4 MeV) is in agreement 



with the experimental data (16.5 ± 0.2 MeV 



I) while the energies of "^Sn (15.1 MeV) and ^"**Pb (12.9 



MeV) are slightly smaller than the experimental values (15.7 ± 0.2 MeV for "^Sn and 13.5 ± 0.2 MeV 



5] 



for ^"^PblSOj). However, the calculated centroid energies for ^^^Sn and ^'"^Pb contain contributions from 
the low-energy discrete states (Fig. which were not considered in the analysis of the experiment. 
Excluding these contributions, both of the centroid energies (15.6 MeV for ^^^Sn and 13.4 MeV for ^''^Pb) 
are also in agreement with the experiments. 

D. Giant quadrupole resonances 



The giant quadrupole resonance is induced by the 2hLL} particle-hole correlation in the harmonic oscil- 
lator potential model, and the isoscalar giant quadrupole resonance (ISGQR) generates a collective state 
in the lower-energy region while the isovector giant quadrupole resonance (IVGQR) generates it in the 
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higher-energy region. In relativistic models the slicU structure is significantly affected by the effective 
mass. As a result, the models with and without the vacuum polarization may yield the different results 
from each other both in the ISGQR and the IVGQR. 

The EWSR distributions of the Coulomb response of the ISGQR (IVGQR) mode are shown in Fig.lTTl 
(Fig. [12) for ^°Zr, "^Sn, and ^"Spb. The ISGQR strength exhibits a strong collectivity. The peaks seen 
in the region of less than 8 MeV for the ISGQR mode come from the low-lying discrete states. In all 
nuclei, the downward shift of the model with vacuum polarization (RHATl) from the model without 
vacuum polarization (TMl and NL3) is clearly observed both in the ISGQR and the IVGQR modes. 
Also, it is found that the RHATl results are in excellent agreement with the experimental data indicated 

3, 3 3 



by the arrow for each nucleus 



We show the RHATl results for the centroid energies of 



the ISGQR in Table IV where they are compared with the corresponding data. The centroid energies 
calculated by the relativistic model with the parameter set NL3 and the relativistic point-coupling 
model were approximately 1-2 MeV above the experimental energies^^ 3- These disagreements are 
drastically improved by the present model including the vacuum polarization. Thus, the effective mass 
plays an important role in the improvement of the results in the ISGQR, while not in the ISGMR with 
compressibility fixed to a certain value. The centroid energies of the ISGQR as a function of the effective 
mass are plotted in Fig. ^1 where we see that the effective mass, as expected, has a substantial effect 
on the energy of the ISGQR. In the present model the large effective mass is theoretically caused by the 
vacuum effect. It is thus shown that the inclusion of the vacuum polarization gives much better results 
for nuclear excitations through the effective mass effect. 

V. SUMMARY 

In summary, we have studied the nuclear excitations by the self-consistent RRPA method with the 
vacuum-polarization contribution extended with the tensor-coupling term for the uj meson. We have 
shown that the vacuum polarization can be easily and duly included in the RPA calculation if the deriva- 
tive expansion for the one-loop nucleon contribution is taken into account in the effective Lagrangian 
and if the vacuum-polarization parts in the unperturbed-polarization tensor are consistent with the 
second-order fluctuation of its Lagrangian. We have indicated numerically that the vacuum-polarization 
contribution is quite important to remove the center-of-mass motion from the isoscalar dipole resonance, 
and the self-consistent calculation involving the Schwinger term conserves the vector current in the 
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multipolc transition. 

Using this method, we have investigated the properties of the coUective excitations of ^°Zr, ^^^Sn, 
and ^°*Pb. With the present model, it is impossible to obtain the effective mass as small as that of 
the conventional relativistic model, since the vacuum polarization works to reduce the scalar-meson 
field. Therefore, the enhancement of the EWSR strength due to the effective mass is suppressed by 
the inclusion of the vacuum polarization. This fact does not contradict with the scattering data and is 
consistent with the analysis of the /3-decay lifetimes. In addition, we stress here that the present extended 
relativistic model is able to give better predictions in almost all giant resonances than the relativistic 
model without the vacuum polarization. In particular, we have shown that the calculated ISGQR 
centroids reproduce the experimental data for each nucleus nicely. Thus the present study suggests 
that the vacuum-polarization effect plays an important role in the excitations through the change of the 
single-particle structure by means of the effective mass effect. Further calculations testing the nuclear 
models for observables with respect to nuclear excitations like a nuclear-polarization correction in muonic 
atoms would be very interesting. 

Finally, we mention the conflict between the relativistic and non-relativistic models for the compression 
modulus. As for the ISGMR in ^"*Pb, the estimation that the compressibility in nuclear matter calculated 
by the relativistic model is in the range K — 250 — 270 MeV is a likely consequence. However, for the 
ISGMR in ^"Zr and ^^^Sn, and for the high-energy part of the ISGDR in the present analysis, the 
compression modulus seems to be a lower value, and therefore the relativistic result may coincide with 
the non-relativistic one. Unfortunately, we could not find a stable solution less than K = 270 MeV with 
the vacuum polarization. In order to resolve this problem, we need the {ui^uj^ f ter m originated by the 
parameter set TMl. The energy of the ISGMR compression modes is not sensitive to the inclusion of 
the vacuum polarization, that is, to the value of the effective mass, which is also the conclusion of the 
non-relativistic study [s^. On the other hand, the energy peaks of the upper component in the ISGDR, 
identified as the compression mode, shift to the lower energy and approach to the corresponding data by 
the inclusion of the vacuum polarization. Therefore, we should study further the relation between the 
compressibility and the effective mass. 
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APPENDIX A: CURRENT CONSERVATION WITH TENSOR-COUPLING INTERACTION 

In order to prove the current conservation for the RRPA transition densities in the present model 
including the tensor-coupling interaction between the u meson and nucleons, we introduce the simplified 
Lagrangian density, 



4m AT 



-[d^A,-d,A,f--(d>^A^)\ (Al) 



where the gauge fixing terms, — 1/2(9^(jj^)^ and — 1/2(9'^ A^J^, are explicitly included. The Noether's 
currents of the photon and to meson fields, 2\ and j^, are same as Eqs. (|30|) and H31() . respectively. Using 
this Lagrangian, the canonical conjugate of the w-meson field reads as, 

rlC f - 

(5^^" - d'^uj^ - ^^a^'^i^ - {d-0JuW\ (A2) 
ocjp zm AT 

which depends on the nucleon field through the tensor-coupling term. 

Now, we shall define the propagator connecting between the photon and the uj meson, Dj^^(x — y); 



where Z?^g(a; — y) and D'^^{x — y) are the free propagator of the photon and uj meson, respectively. Then, 
the four-divergence of the photon-tj meson propagator gives 

U^iDZi^ ~y) = □^(0|T[A^(a;K(y)]|0) === (0|T[4(x)^''(y)]|0), (A4) 

where we use the Euler-Lagrange equation □yl^(.T) = j^(a;). Substituting Eq. I)A3|I into Eq. ljA4|l and 
multiplying -I- m^, we obtain the following relation; 

where we also use the equation (□ -I- rn^)uj^^{x) — jl!j{x). From Eq. ljA2p and the quantization condition 
imposing [j^(a;), 7r^(y)]o — 0, we can easily deduce the relation 

^^^(Q7^7'^ + = {mfA{^)3l{y)m + ^(Olbl(^),^(2/)a'"^V'(2/)]o|0), (A6) 

which is exactly same as Eq. ()29|l . Thus, the straightforward treatment of the Lagrangian including 
the tensor-coupling term naturally introduce the Schwinger term, which is the second term in the right- 

n 

hand side of above equation The Schwinger term plays a crucial role for the current conservation. 
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Multiplying Eq. HA6|I by 9^, one obtains 

(A7) 

The current operators do not commute each other due to the anomalous part, and therefore this 
contribution results in the violation of the conservation law when the conventional polarization tensor 
as indicated in Eq. H17|l only is considered. It can be easily seen that the Schwinger term cancels the 
non- vanishing term of the showing the current conservation of the polarization tensor. 
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TABLE I. Parameter sets of the effective Lagrangians 





RHATl 


RHAT2 


RHAT3 


TMl 


NL3 


m7v[McV] 


938.0 


938.0 


938.0 


938.0 


939.0 


TO<^[MeV] 


390.0 


400.0 


450.0 


511.198 


508.194 


m^[MeV] 


783.0 


783.0 


811.0 


783.0 


782.501 


mp[MeV] 


763.0 


763.0 


763.0 


763.0 


763.0 


9a 


6.046 


6.226 


7.085 


10.0289 


10.217 




8.264 


8.407 


9.220 


12.6139 


12.868 


9p 


10.285 


10.055 


9.155 


9.2644 


8.948 




19.75 


19.39 


18.54 


14.465 


20.862 




-19.42 


-19.42 


-20.82 


3.7098 


-173.31 




0.0 


0.0 


0.0 


71.3075 


0.0 




2.9 


2.8 


2.5 


0.0 


0.0 


A:i.[fm-i] 


1.3021 


1.2993 


1.2894 


1.2905 


1.2995 


^/A[MeV] 


-16.522 


-16.487 


-16.204 


-16.265 


-16.243 


m*/mjv 


0.795 


0.791 


0.774 


0.634 


0.595 


is:[MeV] 


270.8 


281.4 


325.5 


281.1 


271.6 


a4[MeV] 


39.97 


38.68 


34.05 


37.26 


37.41 



TABLE II. Binding energy and charge radius in spherical nuclei 

RHATl RHAT2 RHAT3 TM° NL3 Exp. 

160 

EIA[MeY] "7.976 -7.983 -8.020 -8.043 -8.052 -7.98 

rcft[fm] 2.700 2.699 2.648 2.718 2.686 2.74 

40Ca 

E/A\M.eV] -8.561 -8.557 -8.517 -8.503 -8.547 -8.55 

rchffm] 3.431 3.434 3.398 3.513 3.437 3.45 

90 2r 

E/A[MeY] -8.701 -8.702 -8.656 -8.706 -8.687 -8.71 

rcfe[fm] 4.234 4.243 4.221 4.257 4.241 4.26 

"6Sn 

E/A[Me\] -8.518 -8.521 -8.474 -8.507 -8.485 -8.52 

rcft[fm] 4.576 4.563 4.569 4.605 4.709 4.63 

208pb 

E/A[MeY] -7.859 -7.877 -7.870 -7.875 -7.879 -7.87 

rcfc[fm] 5.489 5.498 5.481 5.489 5.494 5.50 

° The TMl parameter set is supposed to be used above A > 40, hence we write in ^^O and ''"Ca the 
results with the TM2 parameter set, which is arranged to describe the small mass nuclei. 
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TABLE III. Energy-weighted sums of B{EX) over positive-energy RRPA states in unit of e^b'^'-MeV. 
The classical EWSR values (non-rclativistic model) [41] arc also shown for comparison. 





RHATl 


RHAT2 


RHAT3 


TMl 


NL3 


Classical" 


9"Zr 














EO'' 


0.467 


0.471 


0.479 


0.491 


0.494 


0.470 


EV 


3.572 


3.587 


3.655 


4.368 


4.224 


3.303 


E2 


6.843 


6.880 


6.917 


7.315 


7.344 


5.878 


"6Sn 


EO'' 


0.819 


0.824 


0.831 


0.900 


0.902 


0.688 


EV 


4.578 


4.599 


4.599 


5.517 


5.334 


4.229 


E2 


10.59 


10.64 


10.71 


11.32 


11.38 


8.597 


208pb 


EO 


2.177 


2.179 


2.201 


2.374 


2.383 


1.626 


El 


8.044 


8.083 


8.244 


9.680 


9.360 


7.384 


E2 


27.87 


27.95 


27.98 


29.90 


29.98 


20.33 



The radial moments {r'^)p in the classical EWSR are calculated with the charge distribution from the 
MET with parameter set RHAT2. 



TABLE IV. ISGQR centroid-energy positions (MeV). The contribution from the low- lying discrete 
states is eliminated in the RHATl results. 





RHATl 


Experiment 


902r 


14.86 


14.65 ±0.20" 


"6Sn 


13.75 


13.50 ±0.35** 


208pj-, 


11.39 


10.89 ±0.30'' 



" Reference [45] 
'' Reference [51] 
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• Occupied 
o Unoccupied 



FIG. 1: Schematic representation of the polarization tensor: Ho represents a transition from the Fermi sea to 

unoccupied positive- and negative-energy states, while Hf represents a transition from the Dirac sea to unoccupied 
positive-energy states. IIf is divergent without the counterterms. In the present work, its physical contribution 
is evaluated by using the derivative-expansion method. 




2 4 6 8 10 

Excitation energy [MeV] 



FIG. 2; ISGDR strength with the parameter set RHATl in the momentum transfer g = 10 MeV. The spurious 
state gets shifted to zero energy in the full calculation (the solid line). If 11^°^ is neglected in the RRPA calculation 
(the dashed line), the elimination of the spurious state from the real excitation-energy region becomes incomplete. 




FIG. 3: Violation of the current conservation in the isoscalar quadrupole 3.7-MeV state of Pb with the 
parameter set RHATl. The violation is quite large if the Schwinger term is neglected in the polarization function 
(the solid line). The inclusion of the Schwinger term makes the transition current conserved (the dashed line). 
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FIG. 4: The transition-chaxge and -current densities in the isovector antinucleon dipole 1507.2-MeV state of Zr 
with the parameter set RHATl. The transition-current density jio{q) (the dashed line) is considerably large in 
comparison with the the transition-charge density pi{q) (the solid line) and the transition-current density ji2{q) 
(the dash-dotted line). 








1 


1 1 










El transition • 


Zr 








■ 


Sn 








A 


Pb 








~ ~ ~ -* 




































1 


1 







0.7 0.8 0.9 
Effective mass m /mj. 



1.0 



0.7 0.8 0.9 
Effective mass m /m.. 



1.1 



FIG. 5: Energy-weighted sum rules of EO, El, and E2 transitions in ^'^Zr, ^^®Sn, and ^"'^Pb as a function of 
the effective mass m*/mjv obtained with the parameter sets RHATl, RHAT2, RHAT3, TMl, and NL3. The 
classical values axe also plotted in m* /tun = 1- 
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FIG. 6: Calculations of the isoscalar-monopole strength in ^"Zr, ^^''Sn, and ^"^Pb. Arrows indicate experimental 
peak energies . 
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FIG. 7: The centroid energy of the breathing mode in Zr, Sn, and Pb as a function of the compression 
modulus of nuclear matter obtained with the parameter sets RHATl, RHAT2, RHAT3, TMl, and NL3. The 
corresponding data are displayed by the dashed lines |44l |45|. 
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FIG. 8: Calculations of the isovector-monopole strength in Zr, Sn, and Pb. Arrows indicate experimental 
energies p:6. .47^1 . 
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FIG. 10: Calculations of isovector-dipole strength in ^"Zr, ^'^^Sn, and ^"^Pb. Arrows indicate experimental 
energies Isoll . 
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FIG. 11: Calculations of isoscalar-quadrupole strength in Zr, ^^''Sn, and Pb. Arrows indicate experimental 
energies 
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FIG. 12: Calculations of isovector-quadrupole strength in ^"Zr, ^^''Sn, and ^"**Pb. Arrows indicate experimental 
energies |H3iE3l- 
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FIG. 13: Excitation energy of the isoscalar-quadrupole resonance of ^°Zr, "''Sn, and ^°^Ph is plotted as a 
function of the effective mass m*/mjv obtained with the parameter sets RHATl, RHAT2, RHAT3, TMl, and 
NL3. The corresponding data are displayed by the dashed lines [isllSlj . 



